Andreev reflection eigenvalue density in mesoscopic conductors. 



m 

O 

o 

(N 
o 

Q 



p. Samuelsson", W. Belzig^ and Yu. V. Nazarov'^ 
Departement de Physique Theorique, Universite de Geneve, CH-1211 Geneve 4, Switzerland. 
Department of Physics and Astronomy, Universitdt Basel, Klingelbergstr. 82, 4056 Basel, Switzerland. 
Department of Nanoscience, TU Delft, Lorentzweg 1, 2628 GJ Delft, The Netherlands. 

The energy-dependent Andreev reflection eigenvalues determine the transport properties of 
normal-superconducting systems. We evaluate the eigenvalue density to get an insight into for- 
mation of resonant electron-hole transport channels. The circuit-theory-like method developed can 
be applied to any generic mesoscopic conductor or combinations thereof. We present the results for 
experimentally relevant cases of a difi^usive wire and a double tunnel junction. 
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Scattering theory provides a intuitive physical pic- 
ture of transport in mesoscopic conductors P, S 0- 
A large number of transport properties are character- 
ized by the transmission eigenvalues {T^} of the con- 
ductor, the eigenvalues of the matrix product tt\ where 
t is the matrix of the transmission amplitudes and n 
the channel index. Examples are the Landauer formula 
for the conductance G ~ X)™^"' ^^e shot noise power 
P ~ Xln^"'(-'- ~ '^") more general, the full count- 
ing statistics ~ Snl'^ll + rn(exp(ix) — 1)] Knowl- 
edge of the density of the transmission eigenvalues for a 
given mesoscopic conductor allows one to calculate the 
ensemble averaged transport properties. The transmis- 
sion eigenvalue density has been derived for generic con- 
ductors such as diffusive wires 0, Q| , chaotic cavities^ 
double barrier junctions M0|j dirty tunnel barriers 
or combinations thereof [9j . 

When a mesoscopic conductor is connected to a super- 
conductor, the electronic properties of the conductor are 
modified by the induced proximity effect |T2i |. On a mi- 
croscopic level, the proximity effect results from Andreev 
reflection at the normal-superconducting (NS) interface. 
In analogy to the normal conductor, the transport prop- 
erties in the most interesting regime, at energies well be- 
low the superconducting gap, are naturally expressed in 
terms of the Andreev reflection eigenvalues {Rn\- The 
{Rn\ are the eigenvalues of the matrix product Sg^Sg^, 
where Seh is the matrix of the Andreev reflection ampli- 
tudes. The conductance 0| Gns ^ J2n^n and shot- 
noise power Q P/vs ''-^ J2n^n{^ — Rn) are the most 
investigated quantities. 

Due to the dephasing of electrons and holes with an 
energy difference 2E, the Andreev reflection eigenvalues 
depend on energy on the scale of the Thouless energy 
Ec. Only at low energies E <^ Er are the Andreev reflec- 
tion eigenvalues simply related [l5j to the transmission 
eigenvalues as i?„ = T^/(2 — T„)^, providing a direct re- 
lation between the density of Andreev reflection eigenval- 
ues (ARED) and the density of transmission eigenvalues. 
To extend this relation to finite energies, one would need 
not only the correlations between transmission eigenval- 
ues at different energies, but also the energy dependent 
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FIG. 1: Left: Schematic picture of a mesoscopic conductor 
connected to a normal (N) and a superconducting (S) reser- 
voir. Right: The circuit theory representation of two meso- 
scopic conductors: a) a diffusive wire in series with a tunnel 
barrier, b) two tunnel barriers in series. The conductances of 
the respective circuit elements are shown. 



correlations between the transmission eigenvectors, i.e. 
the statistics of the whole energy dependent scattering 
matrix. Apart from being a complicated approach, this 
statistics is known only for a few generic conductors 0|. 
This makes it highly desirable to find a way to directly 
calculate the ARED, without going via the energy depen- 
dent statistics of the scattering properties of the normal 
conductor. 

In this paper we present a theoretical approach to the 
energy dependent ARED for mesoscopic conductors. It 
comes in the form of a circuit theory for quasiclassical 
Green's functions The theory is applied to 

a diffusive wire junction, both with and without a tun- 
nel barrier at the NS-interface, and to a double tunnel- 
barrier junction. In all cases, the induced proximity effect 
leads to an opening of electron-hole channels, i?„ « 1, at 
energies of the order of Ec . Our results provide a unified 
and physically intuitive explanation for a large number of 
theoretical and experimental results on current and noise 
in NS-systems. 

We consider a generic mesoscopic conductor with a 
large number of transverse channels TV ^ I contacted 
to one normal and one superconducting reservoir (see 
Fig. An electron, incident from the normal reservoir 
in channel n, has the amplitude {seh)nm to be backre- 
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fleeted as a hole in channel m. To derive the ARED we 
follow Ref. and define a function in terms of the 

Andreev scattering amplitudes (depending on energy) as 



H{E, C) = Ti' 



ERn 
I - r2 



(1) 



which gives all moments Tr 



when expanded in 



C- From this function, the ARED p{R, E) = (X;„ 5{R - 
Rn)) is given as 



p{R,E) 



ttR^ 



(Im 



H\E,C = 



1 



y/R + iO 



(2) 



where (...) denotes average over impurity (or sam- 
ple) configurations. Next, the Andreev reflection 
amplitudes are related to the retarded anomalous 
Gorkov Green's function F'l'^(Xjjc'), constructed from 
the scattering state solutions 17] to the Bogoliubov- 
de Gennes equation. Introducing F^!^(x,x') = 
Jdydzdy'dz'eniy,z)en{y',z')F^"i^,jc'), with 0niy,z) 
the transverse wavefunction, we get 



(Se/i)„,„ = -ihy/VnVmFl^{x,x') 



(3) 



with F^^{x,x') evaluated at the normal reservoir (x — 
x' = 0) and where v„ is the velocity in channel n. Using 
that F^{x',x) = -[F'f^lx.x')]* , the trace becomes 



Tr 



= -n2^z;„i;„,F^^(x,a;')F„t,(x',x). (4) 



Introducing the generalized quasiparticle current 
/(x,x') = Y.n^'*^riOn{y-, z)en{y' , z')5{x)5{x'), goiug from 
summation over modes to integration over transverse 
coordinates, we can write 



Tr 



(ixi(ix2(ix3(ix4/(xi, X2) 

X i^«t(x2,X3)/(x3,X4)^^^(x4,Xi). (5) 



As pointed out in Ref. the right side of the expres- 
sion can be evaluated at an arbitrary cross-section of the 
normal conductor. To obtain the traces of Sg^jS^^ to all 
orders in Eq. J^l, we introduce a 4x4-matrix Green's 
function in a fictitious ^-dependent field, coupling the 
advanced and retarded components as 



G(X1,X2) = G°(Xi,X2) 4- j dX3dX4G"(xi,X3) 
X C^(X3,X4) [fi + f2] G(X4,X2). 

Here the unperturbed Green's function is 



(6) 



G" = 



G« 
G^ 



qR/A pR/A 



(7) 



suppressing coordinate dependence, and fj are matri- 
ces with the only nonzero element being (fi)42 = 1 and 
(7^2)13 — —1- With these definitions, it follows from the 
perturbation expansion in C, that 



j dxdx7(x,x')tr [fiG(x,x')] ^(H{EX), 



(8) 



where the trace from now on is in the 4x4 matrix 
space. Starting from this expression, we can now ap- 
ply the standard techniques for impurity averaged quasi- 
classical Green's functions, elaborated in detail for the 
circuit theory in Refs. |Sji1^.1S.J- This gives the relation 



tr [fJ{E,C)] /Gq = -iC{H{EX)) 



(9) 



where i{E, C,) is the energy dependent matrix current [l^ 
through the junction and Gg — 2e^ jh. As described in 
detail in Ref. 0|, the terms in Eq. © containing can 
be gauged away from the interior of the junction, giving 
rise to renormalized boundary conditions for the normal 
reservoir 
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N 



(10) 



where G 



RjA 
NO 



±az with (Tz denoting a Pauli-matrix in 



Nambu space. The boundary conditions in the super- 
conductor are not modified, i.e we have Gg^^ — a^- We 
can then, within the circuit theory for matrix Green's 
functions Ref. calculate the matrix current i{EX)- 
Knowing I{E,Q, the ARED can be obtained from Eq. 
(0) and © as 



p{R,E) 



Gq 7ri?3/2 



Re <^ tr 



fil{E,C 



1 



y/R + Tfj 



(11) 

This relation is the main technical result of the paper. 

To illustrate the approach, we first consider a diffusive 
wire with conductance Gat and length L in good contact 
with a normal and a superconducting reservoir. This 
system is described by the quasiclassical Usadel equation 
|23j DdxGdxG = — «£^[(Tz,G], where D is the diffusion 
coefficient. At low energies E <^ Ec — fiD/L? , where the 
proximity effect is fully developed, the density is obtained 
directly from the transmission eigenvalue density of the 
wire giving 
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Pdiff 
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4Gq Ry/1~~R 



(12) 



This is just the same bimodal density as for the trans- 
mission eigenvalues 0, with a prefactor 1/2. In the op- 
posite, incoherent limit, E ^ Ec, the proximity effect 
is suppressed. It was recently shown [l9| that in this 
limit, not only the current and the noise [l^, but all 
moments of the current, can be found by mapping the 
wire onto an effective normal junction consisting of two 
diffusive wires in series. The same mapping procedure 
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FIG. 2: Energy-dependent Andreev reflection eigenvalue den- 
sities (ARED) of diffusive contacts, a) ARED for a diffusive 
wire ideally connected to a normal and a superconducting 
reservoir. Energies are measured in units of the Thouless en- 
ergy Ec ~ fiDjl? and the ARED is normalized to the curve 
for i? = 0. At energies around E^, resonant electron-holes 
channels, i? « 1, open up. b) ARED with a contact resis- 
tance Gt = Gn at the NS-interface. For finite energies the 
ARED has an upper bound; there is a gap in the density for 
eigenvalues above a maximum eigenvalue Rmax- c) ARED 
with a contact resistance Gt = 2Gjv at the NS-interface. For 
E < 10-Bc all eigenvalues contribute, whereas a gap opens up 
for higher energies. 



can be used for the matrix Green's functions giving the 
ARED. As a result the ARED is just given by the trans- 
mission eigenvalue density of the effective junction (with 
conductance G^v/S and the replacement T — > _R, yielding 
p{R, E ^ Ec) — PdiS, which is the same density as for 
low energies. 

For energies of the order of Ec the density differs from 
PdiS- In this regime, the ARED was calculated numeri- 
cally, following the prescription given in Ref. The 
result is shown in Fig. |21 (normalized to the ARED at 
E — Q). Increasing the energy from zero, the density 
of open channels, i? « 1, is enhanced in comparison to 
the density of closed channels, i? « 0, until E k, 5Ec- 
For even higher energies, the density becomes qualita- 
tively similar to pdis again, but with an overall enhance- 
ment [5^ . The induced proximity effect thus results in an 
opening of electron-hole channels at energies of the order 
of Ec- This provides a complementary picture to the ex- 
planation [13 that an interplay of suppressed density of 
states and long-range pair-correlations leads to a mod- 
ification of the transport properties of energies around 
Ec. 

In particular, the opening of additional electron-hole 
channels provides a simple explanation to the reentrant 
conductance peak Moreover, the qualitative differ- 
ence between the regimes E ^ 5Ec and E ^ 5Ec ex- 
plains the suppression of the effective charge qcS^E) = 
(4/3]e(X;„ -Rn(l - Rn))/{J2n ■^») prsdicted theoretically 
in |2l| and found experimentally in 25*1 . The overall en- 
hancement of the ARED for E > 5Ec cancels from the 



effective charge and explains why Qcb = 2e for E > 5Ec- 
In the regime of strong proximity E ^ 5Ec the enhance- 
ment of the density of open channels results, due to the 
factor I — Rn, in a lowering of Coff- 
in the presence of a tunnel barrier (/) at the NS- 
interface (see Fig.^), it was found experimentally that 
due to the proximity effect the zero-voltage conductance 
is strongly enhanced jl^ compared to what was expected 
from classical addition of resistances of the wire and the 
NIS-interface. This was explained in Ref. as reflec- 
tionless tunneling. Backscattering by disorder in the wire 
enhances the chances for Andreev reflection and thus in- 
creases the conductance. Quantitatively, the ARED for a 
tunnel barrier with conductance Gt at the NS-interface 
can be calculated numerically in the same way as without 
the barrier. In Fig. ^jp we present the ARED for different 
energies in the case Gt — Gn- For low energies E ^ Ec 
there is a finite fraction of completely open channels, 
R = 1. Increasing the energy, the whole density shifts 
towards lower eigenvalues and the open channels start to 
close, i.e. there is a gap in the density above a maximum 
reflection eigenvalue i?max < 1 . A similar behavior is seen 
for Gn < Gt, except that the opening of the gap occurs 
for a finite energy. In the example, shown in Fig. [2t for 
Gt = 2Gn, the gap opens for E > lOEc- Note, that 
in the limit Gn ^ Gt where the tunnel barrier can be 
neglected and there is a fraction of open channels for 
all energies. For a dominating tunnel barrier resistance, 
Gt < Gn, there are no completely open channels even 
ior E ~ 0. In the limit Gt Gn the shot-noise was 
recently observed experimentally 28] and studied theo- 
retically [23. A doubling of the full Schottky noise was 
found which can be understood by our results as a con- 
sequence of a concentration of the ARED at small R for 
Gt ^ Gn- 

As a second example, we investigate the ARED in 
a double tunnel-barrier junction (see Fig. Qd), a dif- 
fusive conductor with negligible resistance connected to 
one superconducting and one normal reservoir via tun- 
nel barriers with conductance Gi (to N) and G2 (to S). 
The relevant energy scale for the proximity effect is now 
Ec — {Gi + G2)S/{AGq), the inverse escape time from 
the diffusive region. Here S is the level spacing in the nor- 
mal region. For this system the ARED can be obtained 
analytically. The solution for the Green's function of 
the diffusive conductor for arbitrary energies is given in 
Refs. 30]. Eq. ((TT|l gives after some algebra the ARED 
P = piR,E) 



GlGl 



2nGQR^/'^ 



< R< Ri 
Ri <R< Rii 



(13) 



where p ^ (iGiGi - R(3^) 



■2U/2 



(i?(a2 + /32) - 



AGlGlf/^ a = Gl + Gl- [E/EcfiGi + Gs)^ and 
(3 — 2Gi(Gi -I- G2)E/Ec- The Hmiting eigenvalues are 
Ri = AGlGl/{a^ + P^) and Rn = AGlGj/l3^ (see Fig. 
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FIG. 3: Left panel: ARED for a double-barrier junction with 
G2 = 5Gi for different energies E/Ec- Right panel: Limiting 
reflection eigenvalues 7?/ (solid) and Rii > Ri (dashed) as a 
function of energy for different conductance ratios G1/G2 



13), noting that the region Rj < R < Rjj in the density 
in Eq. H13|) is present only for a < 0. 

The ARED is plotted in Fig. El for G2 = 5Gi, where 
the various properties of the ARED are clearly seen. 
For low energies E < E^[Gl + 01)^1^ /[d + G2), i.e. 
a > 0, the ARED is limited by an upper bound Rj and 
is bimodal, see the upper term Eq. \V6\ . In the limit 
E Ec the ARED reduces to the one obtained from the 
known 9] transmission eigenvalue density for the nor- 
mal double tunnel-barrier junction. For higher energies, 
E > E^{Gl + Glf/'^/{Gi + G2), the lower term in Eq. 

contributes as well. This leads to a qualitatively new 
behavior of the ARED. The peak at the upper bound 
of the ARED splits and an extra bimodal band of An- 
dreev reflection eigenvalues between Ri and Ru emerges. 
Such a peculiar distribution with three singular peaks 
has to our knowledge not been observed before. We at- 
tribute this to a rigidity in position of the electron-hole 
resonaces in double barrier junctions, most pronounced 
for Gi <C G2 . Upon changing the impurity configuration, 
each resonance fluctuates in energy between values re- 
lated to the conductances Gi and G2 . As a consequence, 
for energies below the lowest resonance, a > 0, the reflec- 
tion eigenvalues take on values between and Ri and the 
density is bimodal. For higher energies, a < 0, the peaks 
at and i?/ persist, but there is an additional peak at 
the maximum reflection eigenvalue Ru on resonance. In 
a junction dominated by impurity scattering, this rigidity 
eff'ect is washed out. 

There are two physically different regimes for the 
ARED, depending on the relation between the conduc- 
tances Gi and G2. On one hand, for dominating coupling 
of the diffusive conductor to the normal reservoir, Gi > 
G2, the value of i?max = max(i?/, i?//) never reaches 
unity, i.e. we have no completely open channels. This 
situation resembles the NIS-junction investigated above, 
i.e. the tunnel barrier to the normal reservoir plays the 



same role as the diffusive wire. On the other hand, for 
dominating coupling to the superconductor, G2 > Gi, 
there is always an energy, Ec \/ (G2 — Gi)l(Gi + Gi), for 
which i?inax reaches unity. Of particular interest is the 
limit G2 S> Gi. In this limit, a proximity gap of mag- 
nitude i?c is induced in the diffusive conductor. Only at 
energies very close to E^-, which is exactly at the edge of 
the induced gap, there are now open channels. 

In conclusion, we have provided a theory for the An- 
dreev reffection eigenvalue density in mesoscopic normal- 
superconducting junctions. We find quite generally that 
the proximity effect leads to an opening up of resonant 
electron-hole channels at energies of the order of the 
Thouless energy E^. This was exemplified by studying 
a diffusive wire and a double tunnel-barrier junction. 

The work was supported by the Swiss network MaNEP 
(P.S.), the NCCR Nanoscience and the Swiss NSF 
(W. B.). 
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